General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



facility form 602 


PREPRINT 


A TWO SPACECRAFT TEST OF 
A SINGLE SPACECRAFT METHOD 
OF ESTIMATING SHOCK NORMALS 


R. 

P. LEPPING 

N71 23482 

lACCF’^^lf^W KlUAiiocaL 





^ 

(CODE) 

(NAs« V.K UK tMX or ad NUMBER) 



(CATEGORY) 



MARCH 1971 




GODDARD SPACE FLIGHT CENTER - 

GREENBELT, MARYLAND 


X-692-71-138 



A TWO SPACECRAFT TEST OF A SINGLE 
SPACECRAFT METHOD OF ESTIMATING SHOCK NORMALS 

by 

R. P. Lepping 

Laboratory for Extraterrestrial Physics 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 


March 1971 

Submitted to the Proceedings of the Solar Wind Conference . Asilomar, 
Pacific Grove, California, March 21-26, 1971. 







Abstract 


By assuming the validity of a subset of the Rankine-Hugoniot conservation 
relations for interplanetary (IP) shocks in an isotropic medium it has been 
demonstrated, in principle, that improved shock normals can be calculated 
by using a least squares technique on combined magnetic field and plasma data 
from a single spacecraft. The scheme devised by Lepping and Argentiero (1971) 
uses those six conservation relations not involving pressure and temperature. 
This paper deals with a test of the scheme by examining in detail a shock 
across which the magnetic field changed direction by a small amount (wlO°) . 

On January 26, 1968 at about 1430 U.T. this shock was observed by the plasma 
and magnetic field instruments on Explorers 33 and 35. The spacecraft were 
76.6 and 56.9 earth radii (R ) sunward of the earth respectively (and 43.5 R_ 
from each other), and therefore well outside the earth's bow shock region, 
a necessary condition for a valid test. It was assumed that an IP shock's 
surface is locally plane over dimensions of about 100 R . Using this 
assumption and the known geometrical configuration of the positions of the 
spacecraft with respect to the earth at the times of the shock onset the 
orientation of an "observed" normal was ascertained. For comparison least 
squares best-estimate normals were then calculated for each spacecraft using 
three different time intervals of data in each case: 9, 12, and 18 minutes, 
before and after onset. This was repeated using only the magnetic field 
data and the conventional coplanarity theorem for further comparison. For 
the 18 minute data interval it was shown that the best-estimate normals for 
Explorers 33 and 35 agree with each other within less than 3°, and correspond 
to the "observed" normal within its angular uncertainty due to the time 
uncertainty of the earth's sudden commencement. 


Introduction 
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The purpose of this paper is to test a single spacecraft method of 

estimating shock normals by cross checking the results of the method applied 

independently to the data of two interplanetary spacecraft located about 44 

earth radii (R ) from each other. The method, devised by Lepping and 
£ 

Argentiero (1971), uses a six equation subset of the eight equation Rankine- 
Hugoniot conservation relations for interplanetary (IP) shocks in an isotropic 
medium, i.e., those equations not involving temperature or pressure. They 
showed in principle that improved normals can be calculated by employing a 
least squares technique to best fit the combined magnetic field and plasma 
data from a single spacecraft to three equations of the six equation subset, 
after transformation to an arbitrary frame of reference. The remaining 
three equations are used explicitly to obtain the direction of the normal 
and, provided the average pre-shock plasma velocity is sufficiently accurate, 
the speed of the shock. The reasons for ignoring the equations containing 
temperature or pressure are: 

1. The proton data for these parameters usually show the poorest 

approximation to a step function of all the shock parameters, 

2. Use of these parameters would require electron data which is not 

always available, and 

3. Probably most importantly, use of the energy flux equation, which 

does not take into account possible heat flow across the shock 

front (Hundhausen and Montgomery, 1971), is of questionable validity. 
Observations and Discussion 

This paper deals in detail with a shock whose associated magnetic field 
changed direction by a small angle («10°) across the shock transition zone. 


On January 26, 1968 at about 1430 U.T. this shock was observed by Explorers 

33 and 35 with an 88.8 + 3.6 sec. time delay between them (see Figure 1). 

At about 1441 U.T. a sudden commencement was observed on earth. Approximately 

two hours later (1634 U.T.) Pioneer 8, located about 570 R behind the earth 

£ 

near the earth's tail, saw the shock after some deflection of its normal's 
direction. Notice that the IP shock normal was southward by 20° but at the 
Pioneer location it had become northward by »35° (95% certainty error cone 
angle is 17°). Only Explorer 33 was significantly out of the ecliptic plaiie 

A 

and was 26 R below it. The IP shock normal, n, was almost perpendicular 
E 

o ^ 

(«70 ) to the preshock magnetic field direction, Bj^. Therefore, little 

change of direction of the magnetic field would be expected as the shock 

A A 

passed the two spacecraft. The quantities n and are best fit values, 
whose estimates will be discussed below. 

Figure 2 shows superimposed magnetic field data from Explorers 33 and 35 
around the time of the shock. There was essentially no change in 6 across 
the shock surface and only about 10° change in 0. The horizontal lines 
represent the average of the two individual Explorer 33 and 35 best estimate 
values. The length of these lines indicate the 18 minute time Intervals, 
before and after shock onset, that were used in the best fit calculation. 

All six best fit magnetic field parameters seem to have reasonable values 
when compared to straightforward averages allowing for deviations equal to 
the rms deviations for each. Notice the occurrence of a periodic structure 
before and especially after the shock. Behind the shock the oscillations, 
occurring over about 30 minutes or so, are clearly out of phase between the 
Explorer 33 and 35 observations. 

Figure 3 shows the plasma data also superimposed from Explorer 33 and 35 
observations. The horizontal lines in the pre-shock case are simply averages 
of the dual spacecraft data. However, plasma velocity differences are 
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obtained from the best fit scheme and these along with the added pre-shock 
averages yield the post-shock "best-fit" values shown. Again the lines 
represent an 18 minute Interval before and after the shock. Notice that the 
periodic structure after the shock, which was rather clear in the magnetic 
field data, also appears here except the wave-like signature is not now quite 
as well defined. 

Table 1 gives the best estimate values of the IP shock parameters for 
the two spacecraft and average values of these best estimates. The subscripts 
1 and 2 refer to pre- and post-shock, respectively, and the R-T-N coordinate 
system, centered at the spacecraft of interest, refers to the unit vectors: 

A A A 

R, radially away from the sun in the ecliptic plane; T, perpendicular to R 

AAA 

and lying in the ecliptic such that R x T = N is normal to the ecliptic and 
"northward". W s (Wj^, W^, is the plasma bulk velocity difference Vg-Vj. 
The N's are the number densities and n^^, n^ and n^ are the components of the 
shock unit normal. The Alfven mach nimibers for pre-and post-shock were 8.5 
and 5.5, respectively; these compare well with those of previously studied 
IP shocks (Hundhausen, 1970). The best estimate Explorer 33 and 35 normals 
(calculated from the 18 minute interval) differed by less than 3°. The 
associated calculated shock speeds were 507 km/sec and 520 km/sec respectively, 
giving an average value of 513 km/sec. An average pre-shock plasma bulk 
velocity = (426, 17.4, -7.0) km/sec from the data of both spacecraft was 


used. 


-.352 


Parameters 
18 Minute Interval 


!.) Value 
!xp. 35 

Average of B, 
Exps. 33 and 

-0.24 

-0.92 

-3.07 

-3.12 

-3.83 

-3.66 

-3.41 

-4.28 

-6.09 

-6.35 

-8.26 

-7.90 

85.6 

82.1 

-35.0 

-36.0 

-24.2 

-26.3 

4.45 

4.32 

10.52 

10.1 

0.850 

0.838 

-0.416 

-0.428 

-0.324 

-0.338 
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Two Spacecraft Test 

The best fit IP normal was checked for accuracy by comparing Its angular 
displacement from two fixed and Intersecting lines In space. These lines 
were: first, the segment between Explorers 33 and 35 and, second, that between 
33 and the earth; they Intersected at 47°. Each of these angles can be 
calculated In two ways: first, by a straightforward calculation using the 
best estimate normal; this gives ■ the calculated check-angles, and, second, 
by assuming, for dimensions of about 100 R , 

(a) a plane shock front, and 

(b) a constant shock speed (513 km/sec) and constant normal; 

the latter are the observed check-angles. The calculated and observed 
check-angles can then be compared. 

In the case of the Explorer 33-35 line the observed and calculated angles 
were 80.5° and 84.1° respectively, giving less than a 4° difference. In the 
case of the Explorer 33-earth line the observed and calculated angles were 
47° and 44° respectively, giving approximately a 3° difference. The sudden 
commencement (SSC) time at earth was taken to be 1441 U.T. , giving an 11 
minute delay. If 1440 U.T. Is taken as the SSC time, giving a 10 minute 
delay, the angles become 52° and 44° respectively, an «8° difference. 
Assumptions (a) and (b) for the region from Explorer 33 to the earth, even 
over the path from the bow shock encounter to the earth, are justified within 
an angle error of 8° or so, because the shock "effectively spends" only 
about one tenth of Its total 33-earth travel time In this latter region. 

The (957. certainty) error cone angle associated with the normal was 7.6° 
which Is consistent with the check-angles, or Is perhaps somewhat conservative. 


7 


Comparison of Analysis-Intervals 

To obtain some understanding of the importance of using the proper time 
interval around the shock for the shock analysis, other time Intervals, as 
well as the 18 minute Interval, were used. Henceforth, the term "best 
estimate" refers only to a given analysis-interval for a given spacecraft, 
and not necessarily to the final best estimate of the IP normal. Figure 4 
shows, for the three separate input data intervals, estimates of the January 
IP shock normal, as projected on the R-T plane, for both Explorers 33 and 35. 
The results of both the average magnetic field method (coplanarity theorem) 
and the best estimate method (auxiliary use of plasma data) of estimating 
the IP normal are sho%m. The latter are represented by either dashed or 
solid arrows . and the former by dashed or solid lines ; dashed indicate 
Explorer 33 estimates and solid indicate Explorer 35. The following features 
should be pointed out: 

a. There is a large (70°) spread of the average normals but a reasonably 
narrow (14°) spread of the best estimate normals over the three time 
intervals. 

b. Lengthening the time interval of data around the shock for use in 
calculating the normal does not necessarily improve the estimate, 
even within the short range considered here (i.e. up to 18 minutes). 

c. For each given time interval the best estimate normals between 
Explorers 33 and 35 are closer together than the average normals. 

d. The 18 minute interval was clearly the "proper" choice of interval 
giving a few degrees difference between the Explorer 33 and 35 best 
estimate normals. 

Figure 5 corresponds to Figure 4 except now the estimates of the IF 
normal are projected into the R-N plane. All of the above comnents again 
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hold except for statement "c”. (However, In statement "a" the spread angle 
for the average normals becomes 46*^ but that for the best estimate spread 
remains 14°) . 

If one had been satisfied with only 12 minutes ("B"> Interval) of data 
around the shock and had not taken advantage of the available plasma data 
(or did not have such data) one might have been led Into a false sense of 
certainty about the results because of the relatively good agreement between 
the results of the two spacecraft for this time Interval. 

Conclusions 

We have accurately estimated an Interplanetary shock normal and have 
shown Its direction to be significantly different from the R-dlrectlon both 
In Inclination angle 9 and azimuthal angle 0 (0= ~20°, 0 * 153°); the ecllj^lc 
plane projection was approximately along the average magnetic field spiral 
direction. There was no obvious solar flare associated with this shock. 

The shock may or may not have originated at the sun but It probably did 
not start as a spherical front near the sun unless the front was severely 
distorted over 1 A.U, The periodic structure occurring behind the shock as 
seen especially In the magnetic field data Is no doubt, In part, responsible 
for the fact, stated above, that lengthening the analysis Interval does not 
necessarily Improve the estimate of the normal. A proper anailyslS'lnterval 
Is probably one that encompasses, as exactly as possible, two oscillations 
If such quasl-perlodlc structure exists after the shock, or at most should 
be limited to the Interval Just up to the first obvious discontinuity appearing 


after the shock. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


FIGURE CAPTIONS 

The positions of Explorers 33 and 35 and Pioneer 8, at the 
time of the January 26, 1968 shock, shown in the ecliptic 

A 

R-T plane. Also shown are the best fit IP normal, n, and 
pre-shock magnetic field direction, Bj^, as well as a roughly 
estimated normal (using average magnetic fields) at Pioneer 
8. Quantities in parenthesis refer to the direction 
perpendicular to the ecliptic plane in either degrees or 
R The question mark (?) at Pioneer 8 refers to the large 
uncertainty (error cone 17°) of the normal's estimate at 
that location. 

Superimposed magnetic field data for Explorers 33 and 35. 

F is the magnitude, <t) is the azimuthal angle measured 
counterclockwise in the ecliptic plane from 0 = 0° in the 
direction of the sun, and 0 is the angle of inclination 
measured positive northward from the ecliptic. 

Superimposed plasma data for Explorers 33 and 35. W is the 
thermal speed and V is the magnitude of the bulk plasma 
velocity, whose direction is designated by 6 (same as in 
Figure 1) and ((^ = 0° in antisolar direction) . N is 
the plasma monber density. 

Estimates of the January 26, 1968 shock normal from Explorers 
33 and 35 data and projected into the R-T plane. Both 
average-field and best-fit methods are shown, each for three 


Figure 5 


separate data intervals. 

Estimates of the January 26, 1968 shock normal from Explorers 
33 and 35 data and projected into the R-N plane. 
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